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Abstract. We prove global regularity of solutions of Oldroyd-B equations in 2 spatial dimen- 
sions with spatial diffusion of the polymeric stresses. 



1. Introduction 

We discuss a simple model of polymeric fluid, a complex fluid that is comprised of a solvent, 
which is an incompressible Newtonian fluid, and a dilute suspension of polymeric matter in it. 
The complex fluid occupies a region in physical space M'*, and in this paper d — 2. In models 
f|15). [29]), complicated objects are described by retaining very few degrees of freedom, and the 
polymers are represented by end-to-end vectors m £ M — W^. In general, the configuration space 
M could be much more complicated. We consider as starting point the kinetic descriptions of the 
particles. The state of the particles is determined by a measure f(m)dm, i.e. by a measure which 
is absolutely continuous with respect to the measure dm. The equilibrium measure is obtained by 
minimizing a modified free energy 

(1) £[/] - / (log f + Ui{m) + lu2[f]{m)) fdm 



M 



2 



In this description / log / represents a thermal effect, Ui{m) is a resident potential that describes 
constraints and 

(2) U2[f]{m) = f k{m,p)f{p)dp 

J M 

is a screening (or excluded volume) potential. If the kernel k{m,p) is symmetric and Lipschitz 
then the resulting extremal equation (the Onsager equation), 

(3) / = Z-ie-^[/] 

has always solutions, at least if M is compact ([l3],[14]). Here 

(4) [/[/] =C/i + C/2[/] 

and Z > is a normalizing constant. The modified free energy is not convex, and the Onsager 
equation has rigorously proven phase transitions (when temperature is lowered) and multiple so- 
lutions coexist. The limit of high intensity, (or low temperature) is determined by a popularity 
contest ([H]): a selection mechanism operates by which states which can be easily altered are 
rejected in favor of states whose character is very much like that of their neighbors. The relaxation 
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mechanisms by which the minimum solutions of Onsager equations are approached is given by the 
kinetic equation 

(5) dtf = ediv„i ^/V,„ 

with e a positive constant quantifying inter-particle diffusivity. The equation can be written as 

(6) dtf = eA,„/ + ediv„, (/V,„ ({/[/])) 

and is a nonlinear Fokker-Planck equation. The presence of the nonlinearity shows that the problem 
of deriving this equation from some underlying stochastic system is not trivial. The equation ([SJ 
has £" as a Lyapunov functional. The time derivative of this Lyapunov functional is non-positive, 
vanishes at solutions of the Onsager equation, and if M is connected, only there: 

where 

(8) V^e [ /|V„, (log/ + [/[/])!' dm 

The presence of fluid introduces new degrees of freedom, due to the ambient space, a subset of 
jgid -jpd^ rpj^g particle density acquires x-dependence: / = f{x,m,t). Now the kinetic evolution 
depends on the fluid's velocity u{x,t). This velocity is a function of the macroscopic variables x,t 
alone. The equation ceases to be a gradient equation. The fluid velocity introduces a drift both in 
space and in particle configuration space. The simplest case in which such an effect is immediately 
observed is the case in which a shear transports an elastic rod. The nonlinear Fokker-Planck 
equation is then 

(9) dtf + u- V,/ + (V,u)m • V„J = eA„,/ + ediv„(/V„C/[/]) 

Note that if u is divergence-free in x then {\/xu)m is divergence-free in m S R'^. The ensemble 
of particles might suffer collectively the effect of being in a strained environment. This can be 
modeled by allowing the resident potential Ui to depend on x, t as well as on m. The resident 
potential does not depend on the state /. The self-interaction potential U2[f] depends on the 
macroscopic variables only because / does, but the kernel k is derived from purely microscopic 
information. For more general configuration spaces we have an added drift 

(10) W{x, m, t) ^ V.:^u{x, t) : c{m) 

that depends linearly on the spatial gradient of velocity, and gives a vector field on M. In compo- 
nents: 

W"'{x,m,t) = {dju'{x,t))4-°'{m) 
The coefficients c(m) are smooth. The nonlinear Fokker-Planck equation is 

(11) [dt + u- V,) / -I- div^iWf) = eAmf + ediv™(/V™C/[/]) 

In this equation the coefficients c represent a "macro-micro" interaction that introduces a particle 
drift from the macroscopic drift. The potential kernel k embodies a "micro-micro" interaction, 
while the configuration space AI and the resident potential Ui represent geometric and kinematic 
constraints. The effect that particles might have on the solvent is a more mysterious matter. At 
this level of description, in order to be self-consistent, this effect can only be obtained by averaging 
out the microscopic variables. The effect is embodied in an added stress matrix a{x,t) and it is a 
"micro-macro" interaction, the macroscopic effect of microscopic insertions. An energetic principle 
([!]) states that this "micro-macro" interaction, is such that the coupled system is dissipative. It 
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turns out that this principle is sufficient to provide formulas for the "micro-macro" interaction even 
in non-dilute cases. In the known examples, this principle leads to familiar rules of determining 
the added polymeric stress from the micro-micro and the macro-micro interactions (|18)). In 
(0) this principle is described for non-dilute situations. In general, the fluid velocity solves the 
incompressible (unit density) Navier-Stokes equation 

(12) dtu + u ■ V xU — v/S.xU + V xP — Kdrvxcr 
with K a positive constant with units of velocity squared [err? j sec?)^ 

(13) \lx-u = ^ 

and > the kinematic viscosity. The dimensionless matrix a — {aij{x,t)) represents the added 
stress tensor. We note that only the non-isotropic part of a enters the equation, because changing 
CT to CT -I- AI results in the same equation, with a modified pressure. The coupled system (jlll [T2|) 
is required to dissipate the sum of the kinetic energy and free energy: 



^\u{x,t)\' + K£[f]{x,t) 



dx+ [iy\Vxu{x,t)\'^ + KV{x,t)] dx<0 
In ([3], [5]) it is shown that this requirement (for all initial data) is satisfied if, and only if 

(15) a{x, t) ^ I [c • V„iJ7 - div„ic] fdm 

J M 

and 

(16) {dt + u-\Ix)Ui<0. 

Recall that c is a matrix of vector fields so c • VmU and divmC are linear operators in W^. In the 
case of the Oldroyd B equation, U2 — Q and 

R is positive, M = R'', and 

(18) W{x,t) =Vxu{x,t)m 

i.e. c{m)l'" — 5fmK The formula (1151) gives 



(19) a{x,t)= / {m(E)V,nU)fdm 

(Note that div„c = da{Sfm^) = 0.) 

In the Oldroyd B case a obeys a closed equation, and this is easily obtained by multiplying ^ 
by (m (g) ra) I E? and integrating. The result is: 

(20) [dt+u- Vx)(T = {Vxu)a + ^(V:^^)'^ - 2ka + 2kpl 
where the damping frequency is 

(21) t^±, 
and 

(22) p{x,t) = / f{x,m,t)dm 
Note that from ([9|) it follows that 

(23) idt + u-Vx)p = 0- 
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The mathematical literature on complex fluids is growing, and we cannot give here a complete 
account. Early work (^^, [21]) established local existence results for Oldroyd-B and FENE type 
equations. The Oldroyd-B equations are exact closures of linear Fokker-Planck equations. The 
resident potential is a harmonic potential that is independent of macroscopic variables and the 
self-interaction kernel k vanishes. The "micro-macro" interaction gives cr as a second moment of 
/ and then a self-contained equation for a follows from the linear Fokker-Planck equation. This 
is the only known instance when the Fokker-Planck equation yields a self-contained equation for 
a. The FENE equation is a model in which the resident potential is infinity at a finite extension 
value, prohibiting the particles from extending beyond it. Global existence of weak solutions in the 
presence of spatial diffusion of the polymers was proved in a sequence of papers, the most recent of 
which is (pj). Global existence of weak solutions via propagation of compactness was proved under 
the corotational assumption ([H], [13) ^'^'^i very recently, for the full FENE model ([IH])- There 
is no such result for the Oldroyd B model. The global existence of smooth solutions for small data 
for Oldroyd B-type models was established in ([H], [13) ^^"^ fo'" FENE in ([27^). Global existence 
of smooth solutions for large data in 2D was established for Smoluchowski equations on compact 
manifolds ([S], [Z], [3, [IH], [U, Wl)- Global regularity for large data in the FENE case, under 
the corotational assumption was proved in (11I],[1B]). An approach based on Lagrangian particle 
dynamics was developed in (|22'). Sufficient conditions for regularity in terms of bounds on the 
added stress tensor were established in ([3], [T7]) and further refined in ( 20;). Numerical evidence 
for singularities was provided in ([32j). The paper ([H]) proved global existence for small data 
with large gradients for Oldroyd-B. A regularization obtained by allowing the "spring constant" in 
the harmonic potential to depend on the local rate of strain of the fiuid was obtained also in (|12j). 
Regularity for diffusive Oldroyd-B equations in 2D for large data were obtained in the creeping 
fiow regime in (ilj. More general models will be presented in ([H])- 

2. Diffusive Oldroyd B: a priori bounds 

We have only limited success in the proof of regularity for large data. The case of ([TT|) coupled 
with the 2D Navier-Stokes equations, when M is a compact Riemannian manifold was addressed 
in (|9|) using Fourier analysis techniques and in (TD], [11]) using physical space techniques. If the 
configuration space of the particles is not compact, then the only results of global regularity for 
all data are for modified equations: the corotational case when 2D Navier-Stokes equations are 
coupled via the antisymmetric part of the gradient to linear FENE equations ('27^), the case of 
Oldroyd B coupled with 2D Navier-Stokes when the resident potential responds to excessive rate 
of strain in the fluid ( [H] ) , and the case in which the polymeric stress is allowed to diffuse in space 
([5], in the creeping flow regime). The proofs of global existence for the compact configuration 
space M are based on the fact that a S L°°{dxdt), which easily follows from definitions. The 2D 
Navier Stokes equations forced by the divergence of a bounded stress have unique weak solutions 
that are very well behaved. In particular, they are Holder continuous after an initial transient 
time. This fact was proved in ([lO]). The proof (|lij) of global regularity of the coupled system 
uses then a local existence theorem to address the initial transient time, and to obtain a global- 
in-time uniform bound for the velocity in a Holder space. Then a bootstrap for higher regularity 
is used to finish the proof. The proof of the analogous result in ([9]) uses ideas from ([3]) to 
derive estimates in which the regularity is deteriorating in time, but in a manner that is controlled 
locally uniformly. Assuming a finite time singularity, we reach a contradiction by starting the 
computation close to the putative singularity and showing that the controlled loss of singularity 
forces the solution to remain smooth beyond the blow-up time. The case of M = R'^ with spatial 
diffusion for cr, for Oldroyd B in a creeping fiow regime (coupling to Stokes) in 2D was presented in 
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([B]), and the coupling to the full Navier-Stokes equation is in this paper. For the FENE equations 
(see [5] for numerical results and set-up) ([13) proves local well-posedness in general, and global 
well posedness for the corotational case. The corotational case is one in which the full gradient 
Vu is replaced by its anti-symmetric part in (jlO[) . In that case, if the resident potential is radial, 
then the correct a given by the relation ([T5|) vanishes. If that is the case there is no added stress 
to the fluid, and coupling with (|19p is not energetically balanced. 

The 2D Oldroyd B system is the equation (gOl) coupled with ([H]), We consider the 

variables 

r a(a;,t) = i((7"(a;,i)-a22(x,t)), 

\ c{x,t) = cr"(x,t) -KCT22(a;,t) = Tr(cr(x,i)) 

The matrix cr is symmetric and positive by construction, in view of (|19p. and is given in terms 
of a, b, c by 



(25) 



§+a b 
b f -a 



The positivity of the matrix is equivalent to the positivity of c (which follows from (I19p ) and the 
positivity of the determinant, i.e. 

(26) £! _ - 6^ > 0. 

4 

We denote by 

(27) X{x, t) = i {diu\x, t) - d2uHx, t)) , 

(28) fi{x, i) = i {diu^ix, t) + d2U^ (x, t)) 
and 

(29) uj{x,t) = diu^{x,t) - d2U^ix,t) 

The functions A and fi represent the rate of strain, uj the vorticity. The equations (j20p can be 
written as the system 

Dta = —Lob + Ac — 2ka 
Dfb — uja + /ic — 2kb, 
Dtc = 4Aa + 4^6 - 2fcc + Akp 
Dtp = 

We used the notation Dt = dt + u ■ Wx ■ Let us multiply the c equation by | , the a equation by 
2a, the b equation by 26 and subtract the last two from the first. We obtain 

(31) Dt J - - = -4fc {^-°^ ^ + ^^^PC 
This cancellation of nonlinearity is not surprising because 

(32) — - ~ 6^ = Det (ct) 

and the determinant is conserved along particle trajectories if fc = 0. It is well-known, and easy to 
see that the regularity of the system is decided by whether or not we can bound c in L°° . 



(30) 
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We prove global existence in the case of physical space diffusion of a. The system we consider 
is thus 

(33) (dt+U- V^)cr = oiy^u) + (Va;U)^cr - 2fc(cr - pi) + KAj;Cr 

with K > a spatial diffusivity of the stress, coupled with (fT2|) . (fT3|) . (|23)) . In terms of a, &, c, 
this is the system 

(dt~Vu- Vj;)a = — aj6 + cA — 2A:a + kA^jO 
(9t + u • Va;)6 = + c/i — 2fc6 + kA^;?) 
(9t + u • Va;)c = 4(Aa + ^6) - 2fcc + KAa;C + Akp 
{dt+u-V^)p = 

coupled with (IT^ . p^ . We assume we are in two spatial dimensions, that the initial data for 
a, 6, c, p are in L^(M^) n W^^^^(M'^). We assume that the initial data for the velocity are divergence- 
free and belong to W'^''^{M.'^). We start with a calculation that shows that c is bounded below in 
time, as long as (a, b, c) solve an equation of the form (l34l) with bounded A and /i. It does not 
matter how the functions u, A, p, and uj are obtained, nor even if they are related among themselves. 
We consider the equation obeyed by 



(34) 



(35) 7(a;,i) = c-2v/a2 + 62 
A brief calculation shows 

(36) Da = -2 (k + ) 7 + 4fcp + kA^c - ,J aA^a + bA^b) 
Now we have 

(37) A. (V^^) «^^° + ^^--^ = („2 + ;,2)-i y-(^5^.^ „ ^5^.^)2 > 



so we deduce 

(38) Da > -2 ffc + 4^^^ ) 7 + '^A,7 + 4fcp. 



Va2 + fe^ 

In view of the maximum principle, and the fact that p> 0, it follows that 7(0;, t) > if the initial 
data obey 7(2;, 0) > 0. If the initial data cto is nonnegative, i.e. it has nonnegative determinant, 
do = DcttJo > and nonegative trace cq > 0, then 70 = 2 -^/oq + 6q + rfg — 2-y/aQ + 6q > 0, and 
then it follows that 



(39) c > 2 Va^Tfe2 

for all time. This implies that c > and that (? > 4(a^ + 5^), so we proved 

Proposition 1. A solution of Ii33]) in 2D, advected by a L^{0,T; W'^'°°{M.'^)) incompressible veloc- 
ity, stays nonnegative if it starts nonnegative. 

Remark. This is not an isolated result. In fact any equation of the type 

(40) Dt(T = Ma + <tM'^ + nAa 

with M traceless and bounded, preserves positivity. The proof above is special to two dimensions, 
but the result is true in any number of dimensions, as is easily proved by time-splitting. Indeed, 
if AC = 0, the equation preserves the determinant, and hence positivity. The heat equation, on 
the other hand, preserves positivity because its solution is obtained by convolution with a positive 
function. 
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(41) 



We are ready now to obtain good global apriori bounds. We consider thus the coupled system 

{dt+u- V^)cr = (Vj;u)cr + cr(Va;u)^ - 2k{a - pi) + kA^ct 
(at+w V^)p = 

Vx • -u = 



in R^. We recall that u has dimensions of velocity cm/ sec, p and a are dimensionless, k has 
dimensions of inverse time sec~^, v and k have dimensions of diffusivity cm? j sec, and K of 
cm} I se(? . The equations are dimensionally correct: the first two and the last one have dimensions 
of inverse time and the third has dimensions of cmj scf? . The bounds below can be checked for 
units; this is a way of double-checking that the bounds make sense, and is why we keep the constant 
physical coefficients fc, k, K , without making nondimensional combinations out of them. 

We start by using the energy balance, taking advantage of the fact that cr is a positive matrix 
in view of Proposition [TJ We take the trace of the first equation, integrate in space and add to 
the third equation multiplied by 2u and integrated. We obtain, after the cancellation of the cubic 
terms 

(42) [ + XTr(cr)] dx+ [ [2v\V :,u\'^ + 2fcifTr(cr)] dx < AkK [ po{x)dx 

where po is the initial density of particles. We used the fact that p is transported by a divergence- 
free vector field, and so all its spatial integrals are conserved. We obtain thus 

(43) sup [\\u\\l2+K\\a\\Li]+2iy [ \\Vxu\\l2dt < WuoWh + K\\ao\\L^ + ikKTWpoh^ = Ro 

0<t<T Jo 

We use the fact that c > 2-\/a^ + 5^ and that J^2 cdx is equivalent to the norm of a. The right 
hand side is bounded uniformly, independently of n. We take now the first equation, multiply 
by a and integrate. We obtain 

(44) ^ll'^lli^ + '^llV.^lli^ + 2fc|la|li2 < 2|la|li.|lV.u|U2 + 2k\\po\\ML^ 
Using the Ladyzhenskaya inequality 

(45) \\a\\l,<C\\a\\L4^,a\\L2 
valid in 2d, we deduce 

(46) |||a||i2 + ^llV.alli^ < C [n-^V^uWl^] Ml^ + k\\po\\h 
and, by Gronwall; 

(47) sup ||a||i2 + « / ||V,a||i2 < Ce [||ao||i2 + fcr||po|li2] = i?i 

0<t<T Jo 

It follows that cr e L°^{dt, L^(R'^)) and a e L'^{dt, W^''^(R'^)) are bounded apriori in terms of the 
initial data. This used the previous bound (j43p . and that used the positivity of a. Taking the curl 
of the third equation in (|41l) we obtain 

(48) dtUJ + u- V^uj - vl^uj = ■ ^\N^a 

Multiplying by w and integrating we obtain, after integration by parts and Young's inequality: 

(49) ^ll"""^^ + ^\\^M?L- < ^l|V.a||i2 
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Using (1471) we have an apriori bound for the time integral of the right-hand side, and so we obtain 



(50) 



suPo<t<T 11^11^2 +1^/0 llV^wll^^ 

II ^2 +-«^ II "0 II i 1 li PO II i 1 



CK 
= i?2 



[Ikol 



2 



fcr|ipolli2] + ||c.olli2 



Now we take the first equation in (|4T)) multiply by — Aa;(T and integrate. We obtain 

(51) 2^ll^=^^lli^ + '^ll^^^^lli^ < ||V^M||i2||V^fT||2, + 2 [||(7||i4||V^u||i4 + k\\p\\L2] ||A:,cr||i2 

The first term in the right-hand side comes from the advective term after integration by parts and 
use of incompressibility. We bound it using Ladyzhenskaya's inequality and Young's inequality. 
We arrive at 

-\W\\l4VM\h ' '^^ ' 



at K 



\Po\\h 



K K 

From Ladyzhenskaya's inequality, (^5]) applied to a and to V^jU, from (|T7)) and ([SUl) and usual 
elliptic estimates relating V^^u in to to in the same space, we see that 

(53) 



W\\l4^xu\\l4 < B 



where the bound B is an explicit expression in terms of the initial data, coefficients and T: 

C 



(54) 



B = 



R1R2 



where i?2 are the right-hand sides of (|T7)). ((50)) . Note that B is dimensionless. It follows then 
that 



(55) 



"0llr2+'^ll''0llrlH-4fc^fT|IP0llrl 



< Ce- 



||V.^olli2+i3 + ^||polli2 



Armed with this, we return to (1481) . multiply by AxU) and integrate in space. We obtain, after 
usual manipulations 



(56) 



|||V.H|i2+H|Ac.||i2<cM^ 



We use the interpolation inequality 

(57) hlli^e <C||u|U2||Au|U2 

valid in and thus, using (j43|) . (|50p and ([55]) we obtain 

supo<t<T l|Vxw|||2 +v ll IIAwlllsdt 



C— IIA.alli. 



(58) 



|V:rWo| 



< Ce^" 'VTfl:o-R.2 
We differentiate ^ 

(59) dtiy^p) ^ u ■ Vxiy^p) = -iy^ufv^p 

We use the interpolation inequality 

(60) W^xuWl^ < C||V,u||i2||AV,u|lL2 



— i?4 
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to deduce that 

(61) sup |1p|1h/i,2(r2) < e'^-*^^'^*^! ||po|ki,2(R2) = i?5 

0<t<T 

We proved: 

Proposition 2. Let u,p, a,p be a strong solution of |^7p on [0,T], with ao G L^(M^) n W^'^{M^^) a 
nonnegative symmetric matrix, po G L^(M^)nVF^^^(R^) nonnegative anduo € M^^'^(M^) divergence- 
free. 

Then p e L°°(0,r,Li(M2) n W^''^{R^)), u e L°°(0,r,VF2'2(M2)) nL2(o^T, W^3.2(k2)) ^ ^ 
L°°(0,T,Li(E2)nW^i'2(i^2))p^2(o^r^^2,2(jfj2>)^ ofcei/ t/ie hounds 0, g^j, fSgj, (igj, {5^, and 

The bounds above used the equations and are fragile: they use the positivity of the matrix a. 
There are several methods of approximating solutions, but these (Galerkin approximation, mild 
formulation and Picard iteration) do not preserve positivity. The simplest proof of global regularity 
is therefore done then by continuation. We show (by any method) that a strong solution exists 
and is unique for a short time. Because it solves the nonlinear equation, if a starts positive, it 
stays positive by Proposition [U and then, by the a priori bounds, the solution has controlled size 
at the end of the short time existence. Existence and uniqueness of strong solutions then allows 
to conclude that the maximal time of existence is infinite. 

3. Diffusive Oldroyd B: Local and Global Regular Solutions 

Proposition 3. Let uq e M^^'^(M^) be a divergence- free vector. Let ctq G W^''^{M.'^) be a symmetric 
matrix, and let po £ L^(M.'^) n W^''^{M.'^). There exists a time Tq > depending on the norms of 
Mo,<TojPo o-i^d a unique strong solution 

U e L°°(0, To, PV72'2(]R2)) n ^2(Q^ PV73'2(K2)) 

(62) a G i°°(0, To, VKi'2(R2)) n ^^(o, To, W^-^{R^)) 

p e l°°(o,To,l1(]r2) n M^i^2(j(g2)) 

of 1-^ j[ ) with initial data UQ,ao, Pq. 

Above P is the Leray-Hodge projector on divergence-free vectors. We sketch a proof via an 
explicit scheme, a Picard iteration for the mild formulation of the equations. We take the Banach 
space defined by the right hand sides of with the corresponding norm || (u, cr, /9)||b (the sum 
of the respective norms). 

\\{u,a,p)\\B = \\u\\x + h\\Y + \\p\\z 
where X = L°°{0,Tq,FW^'^{R^)) n L^{0,To,FW^'^{R^)), 

Y = L°°{0, To, W^-^{R^)) n L^{0, T, W^-^{R^)), 

\W\\y = ||cr||L~(0,To,IVi.2(R2)) + ||o'||l2(o^To,W2.2(r2)) 

and Z = L°°{0,To,L^{R^) f] W^^^{M.^)) with norm 

\\P\\Z = sup (||p(i)||Ll(R2) + ||p(t)||vi/l,2(R2)) 

0<t<To 



10 



PETER CONSTANTIN AND MARKUS KLIEGL 



new „new ^new 



We set up a fixed point equation U = F{U) in B for tfie triplet U = {u,(j,p), where F{U) 
(63) 



) given by 

new 



a 



C-^uq + Qi{u, u) + Li{(j) 



[t) = N{u), where {dt + u ■ V^)p^'''^ = 0, 



p"^-(x,0)=po(a;) 



Jo 

L,{a){t)=K f 
Jo 



P(div^cr(s))ds 



L2{p){t) = 2k e(*-")('''^-2'=V(s)Ms 



where 
(64) 

(65) 

(66) Q2iu,a){t) 
and 
(67) 

We check that 

\\Qi{u,v)\\x<e\\u\\x\\v\\x 
\\Li{a)\\x < C\\a\\Y 

CggN ||Q2(u,cr)||Y < e||u||x||cr||y 

ll-^2(p)||y < e||p||L~(0,To,L2(K2)) 

||^(u)||z < C'||po||Li(R2)nH'i.2(R2)e" 
\\Niu) - N{v)hinL- <4u- vWxe'^^-Mx+WvW^) 

where e can be made as small as we wish by choosing To small enough (while keeping everything 
else fixed) and C is a fixed constant. Obviously Qi and Q2 are bilinear, Li and L2 are linear, 
and the map N is nonlinear. The inequalities (1681) require little effort to check, and they involve 
straightforward inequalities including the Ladyzhenskaya inequality and maximal regularity of the 
heat equation. For instance, the function 

(69) qi{x,t) ^ Qi{u,v){x,t) 
obeys 

(70) 9tgi - :/A,<7i = P(u • V,i;), gi(x,0) = 
and the bound is easily obtained by energy methods, via 

d 



2dt 



\Aq 



i^\\V^Aqi\\i2 < \\V^{u-V^v)\\L2\\Va:Aqi\\L2 



the Ladyzhenskaya inequality for the term involving V^u^xV and the interpolation inequality ([57 
for the term involving u ■ Vx(Vj;w), to yield: 



(71) 



SUPo<f<T 11^91 



< 



^/„^''[||V.^|U.||A«|U.||V.t;|UH|A^;|U^ + ||«IUH|Au|U.||Az.||i.]di 



L2(R2) 



i|li2(R2)rfi 



<CTo\\urx\\vrx 
Next, li{x,t) — Li[a){x,t) obeys 



(72) 
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Applying A and multiplying by A^i we obtain easily the desired estimate 

(73) sup ||A/i||i. / llV^A^illi. < / \\Aa\\l,dt 

a<t<To Jq V 

The estimate for 172 — Q2(w, cr) is straightforward and is easily obtained by multiplying the equation 

(74) dtq2 - KAq2 + 2kq2 = -u ■ V^cr + (V^u)cr + a(y^uf 

by Aq2 and integrating. The estimate for I2 — L2{p) is obtained in the same fashion. The bound 
for N{u) is simply obtained by differentiating. The bound for N{u) — N{v) in D uses the a 
priori bound for individual N(u) and N{v). 

We set up the iteration f/"^^ — F{U""'). By choosing e small enough it is then easy to see that 
the sequence C/" is bounded in B and converges exponentially fast. Let us verify first that we can 
prove by induction that there exist constants A, F, D and e such that 

(75) < r, 

llp"lk<^ 

By induction we have 

\\u"+^\\x<Ao + eA^+Cr, 
<ro + eAr + ei?o, 

where AQ,ro, Rq, Ri depend only of the norms of the initial data. We can choose for instance 
A = 2{Aq + CT) and solve first the induction inequality 

r > To + 2e{Ao + Cr)r + eRo 

by choosing e small enough. Then, by choosing e smaller we guarantee that eA < ^ which implies 
that 

A>Ao + eA'^ + CT 

Then D is chosen in function of A. The induction succeeds because ||p"||l°°(o,To,l2(r2)) < Ro is 
bounded a priori. Once we know that the sequence is bounded in B, then we can obtain exponential 
convergence of u" in X, a" in Y and p" in n i^. 

Theorem 1. Let uq e Vl^^^^(M^) he divergence- free. Let po G L-^(M^) n W^'^(R'^) be nonnegative 
and let ctq G L^{M.'^) f] W^''^{M.'^) be a nonnegative matrix. Let T > be arbitrary. There exists a 
unique strong solution of J/^ j| ) satisfying the bounds |^7| ), i50\) . i55\) . i58\) . and H61\) . 

The proof follows along the lines described above: by Proposition [3] a strong solution exists 
for a short time interval [0,To]. We consider the maximal interval of existence: Ti — supTo < T 
such that the solution is strong on [0,To]. By Proposition [T] the matrix a remains positive. By 
Proposition[2]the solution obeys a priori bounds that are uniform on [0, T]. (The a priori bounds are 
bounded uniformly for Tq <T). Then, it must be that Ti = T, because otherwise, by Proposition 
13] we could extend the solution beyond Ti, contradicting its definition. 

Theorem 2. If uq G W"^+'^''^{R'^) is divergence-free, CTq G L^{R'^) n VF"'2(M2) is a nonnegative 
symmetric matrix and Pq G -L^(]R^) n W'"'^(R^) is nonnegative, m> 1, then the solution obeys 

u G i°°(0, T, IV^+^'^iR^)) n L^{0, T, VF"+2'2(M2))^ 

cr G L°°(o,T,ii(R2) n w"^^^{M.^)) n i^(o,r, w'^+^^^iR^), 

p G L°°{0,T,L\R^) n W'^'^iR^)). 
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Remark. The higher regularity assumption on the initial data is needed because p does not 
obey a diffusive equation and is not constant. If po = 1 then p = 1 for all time. In that case, it 
follows from parabolic regularization that u,a G C°° {to,T,M.'^) for positive to and the initial data 
at i = can be much less regular. 

4. Conclusion 

We proved global regularity for the Oldroyd B system coupled with the full incompressible 
Navier-Stokes system in 2D when the polymeric stress is diffusing in space. This kind of model 
was considered for weak solutions in several studies ([!]) and was used in numerical investigations 
([32], [32] )• Our results use the positivity of the stress matrix, which is preserved under diffusive 
evolution. 

Acknowledgment. PC's research was partially supported by NSF-DMS grant 0804380. 
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